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Jona-Lasinio model with the effective AS = 1 nonleptonic weak interaction. 
The AI = | amplitude is in reasonable agreement with experimental data. 
On the other hand, the calculated AI = | amplitudes strongly depend on the 
mass of the low-lying scalar-isoscalar a meson, and therefore give a strong 
constraint on the parameters of the model. 

I. INTRODUCTION 

The AI = | rule in the K — > 2tt decays is one of the interesting problems of the low- 
energy hadron physics. Since all the particle in the initial and final states in these decay 
processes are the octet pseudoscalar mesons, which are considered as the Nambu- Goldstone 
(NG) bosons associated with the spontaneous breaking of chiral symmetry, one can use all 
the known properties of the NG bosons in the analyses of the processes. 

One of the widely used frameworks to study the low-energy properties of the NG bosons 
is the chiral perturbation theory (ChPT) Q. In ChPT, we introduce an effective lagrangian 
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constructed in terms of nonlinearly transforming NG boson fields and treat the explicit 
breaking of the chiral symmetry and the external momenta of the NG bosons perturbatively. 

In order to study hadronic weak interactions at low energy, an effective weak interaction 
lagrangian is commonly introduced. For K — > 2n decays, we use AS* = 1 nonleptonic weak 
interaction lagrangian Cyy, which has been derived from the standard theory by integrating 
out heavy degrees of freedom. The renormalization group equation with one-loop QCD 
corrections is solved to give Cw at a hadronic mass scale /i in terms of local four quark 
operators 

Using Cw, the K — > 2n decays have been studied in the ChPT J7]||. In this approach, 
the four-quark operators in C w have to be mapped on to the operators of the NG bosons. 
Thus, one should examine whether this treatment is justified. It should be also noted that 
the K — > 27r decay amplitudes are expected to be sensitive to the explicit breaking of the 
chiral symmetry because the matrix elements vanish in the 577(3) limit. In ChPT, this 
symmetry breaking is represented by the finite masses of the NG bosons. It is, however, 
controversial whether the ChPT can qualitatively handle the explicit symmetry breaking in 
the strangeness sector jnj. Therefore it seems interesting to study the K — > 2-k decays in 
an effective theory with the symmetry breaking in a direct form. 

In order to investigate these two points, it may be a good idea to use a chiral effective 



quark model such as the Nambu-Jona-Lasinio model [p[| . If the quark fields in the chiral 
effective quark model are identified with those in C w , there is no ambiguity of using C w . 
Another good point of the chiral effective quark model is that the explicit breaking of the 
chiral symmetry can be introduced naturally by the current quark mass term. 

A few years ago, Morozumi et al. [12| derived a chiral weak lagrangian by performing 



the bosonization of the Cw using the NJL model as a guide and applied it to the K — > 2-k 
decays. They compared their results with those in ChPT and found that the AI = | 
amplitude receives additional enhancement and thus the long-standing AI = | enhancement 
problem is solved. However in ref. the derivative expansion was used and the current 



s-quark mass was treated as a free parameter. Furthermore they calculated only the leading 



contributions in the expansion. 

We present fuller calculations of the K — > 2n decays in a Ul{3) x Ur(3) NJL model 
combined with the effective AS = 1 nonleptonic weak interaction lagrangian Cw in the 
Feynman diagram approach. Here the low-lying meson states are calculated by solving 
the Bethe-Salpeter type equations in the ladder approximation without using the derivative 
expansion. The model parameters, such as the quark masses, are chosen so as to be consistent 
with the observed NG boson properties. We calculate not only the leading terms but also 
the next order terms in the expansion of the K — > 2tt decay amplitudes. 

Although the NJL model does not confine quarks, it is considered that the NG bosons are 
strongly bound and therefore the processes containing only the NG bosons, such as K —>■ 2n 
transitions, can be described by this model fairly well. It is also known that the low energy 
theorems of the NG bosons are satisfied in this model. 



II. EXTENDED NAMBU-JONA-LASINIO MODEL 



We work with the NJL model lagrangian density extended to £/z,(3) x Ur(3) case HlO| : 



Cnjl = Jtidrf - m)il> + E[0M» 2 + • (1) 

1 a=0 



Here the quark field ip is a column vector in color, flavor and Dirac spaces and A a is the 
£7(3) generator in flavor space. The lagrangian Cnjl incorporates the current quark mass 
matrix m = diag(m u , m s ) which breaks the chiral Ul(3) x Ur(3) invariance explicitly. 

Quark condensates and constituent quark masses are self-consistently determined by the 
gap equations. We assume the isospin symmetry, i.e., m u = m^. The pseudoscalar chan- 
nel quark- ant iquark scattering amplitudes are then calculated in the ladder approximation. 
From the pole position of the scattering amplitude, the pion mass m n and the kaon mass 
rriK are determined. The pion-quark coupling constant g n and the kaon-quark coupling con- 
stant gx are determined by the residues of the scattering amplitudes at the pion and kaon 
poles respectively. The pion decay constant /„- and the kaon decay constant fx are also 
determined by calculating the quark-antiquark one-loop graphs. 



III. EFFECTIVE NONLEPTONIC WEAK INTERACTION 



The effective AS = 1 nonleptonic weak interaction lagrangian density we have used is 

= % E (<% + G %) Qr (2) 

with £i = VidV*. and V is the Kobayashi-Maskawa matrix. The four-quark operators Q r are 



Ql = {s a d a )y_ A {U[ 3 Up) v _ A , Q 2 = {s a d f3 )y_ A {llpUa) 



aJy-A ' 

Q3 = {s a d a ) v _ A Y (Qi3Qi3) v _ A , Qs = {s a d a ) v _ A Y {Wp) v +ai 

q=u,d,s q=u,d,s 

(3) 

q=u,d,s 

where a and (3 are indices in the color space. Qi and Q2 have both the AI = 1 and AI = | 
components while Q3, Q$, Qq and Q2 — Qi have only the AI = | component. As for the 
Wilson coefficients C° and C*, we take the results at the energy scale /1 = 0.244GeV with the 
top quark mass m t = 200GeV in ref. |§. The numerical values are C\ = —0.216, C2 = 0.346, 
C 3 = -7.12 x 10~ 3 , C 5 = 3.89 x 10~ 3 and C 6 = -0.021 for C r = C% + C%. 



IV. K -> 2vr DECAY AMPLITUDES 

The invariant amplitude for the K — > 2-n decay Tk^i-k is defined by 

( vr^M^) 1^(9) ) = (2tt) 4 5 4 (A; 1 + k 2 - q) T K ^. (4) 
By calculating the Feynman graphs shown in fig. 1, we obtain the following results. 

T K +^+^ = -J= (Ci + C 2 ) 7y| ( X + ]y~) X ' ( 5 ) 

Gf ( C\ C s C5. 



T^o^o = (-^X - C7 2 X + + j^-Y + C 6 Y ) . (7) 



Here X and Y = (Yi + F 2 + F 3 ) are 
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X = -2y/2&g K 



d 4 p 



T r (c,D) 



d A p 

W) 4 



Tr^ \s s f (p - hh»s u F ( P h 5 s d F (p + fefre 



(8) 
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d 4 p 



Tr^ \S F (p) l5 S d F (p + k 2 ) l5 



d A p 



sUp-k^sUpHsUp + hH 

Y 2 = (l + T a (q 2 )) I ^Z Tr(C ' D) [ S f(p - qHS d F (pH 

_T r (c,D) 



(9) 



d 4 p 



(2tt) 

S d F (p + q)S F (p)^S F (p + k 1 )^ 



(10) 



Y 3 = -Ay/2fig K J ^ {Tt^ [S F (p)} - Tr^ [S F (p)}} 



x 



d A p 

W) 4 



tt<*» \s F ( P ) l5 s F (p) l5 s F (p + h) l5 s d F (p + q ) l5 



(ii) 



where Tr^ c ' D ^ means trace in color and Dirac spaces. The effect of the gg-collective a meson 
in the intermediate state is expressed by 

G S I(q 2 ) 



T a (q 



1 " G s I(q 2 ) 



(12) 



with 



Hq 2 



1 ' ^ T r {c ' f ' D) [SFm u SF(p + q)\ u } 



(13) 



where Tr^^ D ^> means trace in color, flavor and Dirac spaces and A" = y §A° + y |A 8 is a 
matrix in flavor space. Fig. 1 (e) does not contribute because the integral 

d 4 p 



(2vr) 4 

vanishes in the isospin symmetry limit. 



Tr^ S u F (p + q)YS F (p) l5 S d F (p + h) l5 



V. NUMERICAL RESULTS 



The experimental values of the invariant amplitudes for the K — > 2n decays are as follows 
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\T K +^ +7T o\ = 1.83 x l(T 8 GeV, 
|2*WNr-| = 2.76 x 10~ 7 GeV, 

\T K o^w\ = 2.63 x 10~ 7 GeV. (14) 

Here the K + — > decay is a pure AI = | decay and the other two decays have both 

AI = | and AI = | components. 

In our theoretical calculations, the parameters of the NJL model are the current quark 
masses m u = rrid, rn s , the coupling constant G$ and the covariant cutoff A. First, we take 
m u = rrid as a free parameter and study the m u dependence of the K — > 2ir decay amplitudes. 
Other parameters m s , Gs and A are determined so as to reproduce the observed values of m n , 
rriK and f n . The calculated results of \T K +^+^\, \Tk°->-k+k-\ and IT^-o^o^o | are shown in 
fig. 2. As can be seen from fig. 2, the AI = | amplitude has little dependence on m u , while 
AI = | amplitudes strongly depend on m u . This strong dependence is attributed to the T a 
term in Y2, eq. (10). In fig. 3 we show the calculated constituent quark masses M u and M s as 
functions of the current u-quark mass m u . As a result of our fitting procedure, the constituent 
quark masses almost linearly increase when the current u-quark mass increases. The quark- 
antiquark scattering amplitude in the scalar-isoscalar channel has the pole of the a meson 
just above the constituent quark- ant iquark threshold in the NJL model Therefore 



T a (q 2 = m 2 K ) quickly increases when m a (~ 2M U ) becomes close to mx by changing m u . 

If we fit m u so as to reproduce the AI = | amplitudes reasonably, then we obtain the NJL 
model parameters, m u = m d = 7.1MeV, m s = 178.2MeV, A = 0.844GeV and = 0.66. 

The calculated constituent quark masses are M u = Md = 286.2MeV and M s = 529.3MeV 
and the meson-quark coupling constants are g w = 3.023 and g K = 3.254. We have used 
fn-n-, friK and /„- as inputs, so fx is the prediction and our result is fx = lOlMeV, which 
is about 11% smaller than the observed value. The ratio of the current s-quark mass to 
the current u,d-quark mass is m s /m u = 25.1, which agrees well with m s /m = 25 ± 2.5, 
(m = \{m u + md)) derived from ChPT |T5| . 

The calculated results of the K —>■ 2tc decay amplitudes in the above parameter set are 
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|7k+_^o| = 2.52 (1.89) x l(T 8 GeV, 
\T K o^+^-\ = 2.65 (2.95) x 10~ 7 GeV , 

17^0^0^1= 2.30 (2.68) x l(T 7 GeV . (15) 

Here the numbers in parentheses are the results without including -^U terms. The 
terms enhance the \T K +^ n + n o\ and reduce the \T K o^ n + n -\ and \T K o^ n o w o\. This is opposite 
to the nonlinear sigma model case. It is possible that the full QCD contains other 
corrections which are not taken into account in the NJL model. The AI = | amplitude is 
our prediction and it is about 38% bigger than the experimental value. In ref . fl2| , K —>■ 2ir 
decay amplitudes have been calculated in the leading order of the derivative expansion 
without the terms and their result is |'7^+_ j . 7r -i- 7r o | = 2.5 x 10 -8 GeV. So the higher order 
terms of the derivative expansion improve the agreement to the experimental values. 

The contributions from each diagram are as follows. X = 2.49 x 1(T 2 GeV 3 , Y 1 = 
0.274 GeV 3 , Y 2 = -1.564 GeV 3 and Y 3 = -4.1 x 10" 4 GeV 3 . The contribution of Y 3 type is 



not included in the ChPT and the bosonization approach in ref. As far as we know, this 
is the first estimation of this type of the contribution, however our numerical result shows 
that its contribution is negligible. 

The main origin of the enhancement of the AI = ^ amplitudes is the existence of the low- 
lying scalar-isoscalar a meson, which is not experimentally confirmed. It is argued that the 
cr meson has a very large a — ► 2tt decay width and thus is not observed. In our calculation, 
the a — > 27r decay width is not taken into account although the a meson has a small (50MeV 
~ lOOMeV) a — > qq decay width. We estimate the effect of the strong a — > 2tc decay on the 
K — > 27r transitions by using the following scalar-isoscalar form factor 

( 2M„-il) -q 2 



' 2 , 

instead of T a (q 2 ) in eq. (12). Here T is the parameter which represents the a —>■ 2n decay 
width and g a is determined so as to reproduce our results given in eq. (15) in the case of 
r = 0. If one takes T = 500MeV, then the amplitudes become \T K o^+-n- 1 = 8.59 x 10 _8 GeV 



and iT^o^TrOTrO | = 6.52 x 10 _8 GeV, which are about | of the original values. Even in this case, 
if the model parameters are refitted again, one can reproduce the K — > 2ir decay amplitudes 
reasonably by taking m a close to mj-. 



VI. CONCLUDING REMARKS 



Using a three-flavor Nambu-Jona-Lasinio (NJL) model with the low-energy effective 
AS = 1 nonleptonic weak interaction, we have studied the K — > 2tt decays. As suggested 
in ref. |12|], we have found that the existence of the a meson with the mass close to nix 



enhances the AI = | amplitudes and all the K — > 2tt decay amplitudes are reproduced rea- 
sonably well. The role of the a meson in the low-energy hadron physics is not established so 
well. In the K — > 2ir decays, the weak interaction gives rise to the flavor change s — > d and 
it is one of the reason why the a plays an important role in the processes. It may be useful 
to study the K — > 2ir decays and the tctt scattering processes simultaneously by including 
the two pion intermediate states in the analysis since the a meson plays an important role 
in the tctt scattering processes too. It is left as a further study. 
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Figure Captions 



Fig. 1 Feynman diagrams for the K — > 2ir decays. The dotted lines correspond to the 
gg-collective mesons in the NJL model while the solid lines to the constituent quark 
propagators. The hatched ellipses represent the effective AS = 1 nonleptonic weak 
interactions. The diagrams (a), (b), (c) and (d) correspond to the X, Y 1: Y 2 and Y 3 
contributions respectively. 

Fig. 2 Dependence of the calculated K — > 2n decay amplitudes on the current u-quark 
mass. The solid line represents 10 x T K +^ n + n o and the dashed line corresponds to 
T k o^+-k- while the dash-dotted line to Tr-o^o^o . 

Fig. 3 Dependence of the calculated constituent u-quark mass M u and the constituent 
s-quark mass M s on the current u-quark mass m u . 
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This figure "figl-l.png" is available in "png" format from: 
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Fig. 1 

(a) 
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Fig. 2 
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This figure "figl-3.png" is available in "png" format from: 
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Fig. 3 
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